A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

The association of obesity-related traits on COVID-19 severity and hospitalization is affected by socio-economic status: a multivariable Mendelian randomization study.

BACKGROUND: Due to its large impact on human health, socio-economic status (SES) could at least partially influence the established association between obesity and coronavirus disease 2019 (COVID-19) severity. To estimate the independent effect of body size and SES on the clinical manifestations of COVID-19, we conducted a Mendelian randomization (MR) study. METHODS: Applying two-sample MR approaches, we evaluated the effects of body mass index (BMI, n = 322 154), waist circumference (WC, n = 234 069), hip circumference (n = 213 019) and waist-hip ratio (n = 210 088) with respect to three COVID-19 outcomes: severe respiratory COVID-19 (cases = 8779, controls = 1 000 875), hospitalized COVID-19 (cases = 17 992, controls = 1 810 493) and COVID-19 infection (cases = 87 870, controls = 2 210 804). Applying a multivariable MR (MVMR) approach, we estimated the effect of these anthropometric traits on COVID-19 outcomes accounting for the effect of SES assessed as household income (n = 286 301). RESULTS: BMI and WC were associated with severe respiratory COVID-19 [BMI: odds ratio (OR) = 1.51, CI = 1.24-1.84, P = 3.01e-05; WC: OR = 1.48, 95% CI = 1.15-1.91, P = 0.0019] and hospitalized COVID-19 (BMI: OR = 1.50, 95% CI = 1.32-1.72, P = 8.83e-10; WC: OR = 1.41, 95% CI = 1.20-1.67, P = 3.72e-05). Conversely, income was associated with lower odds of severe respiratory (OR = 0.70, 95% CI = 0.53-0.93, P = 0.015) and hospitalized COVID-19 (OR = 0.78, 95% CI = 0.66-0.92, P = 0.003). MVMR analyses showed that the effect of these obesity-related traits on increasing the odds of COVID-19 negative outcomes becomes null when accounting for income. Conversely, the association of income with lower odds of COVID-19 negative outcomes is not affected when including the anthropometric traits in the multivariable model. CONCLUSION: Our findings indicate that SES contributes to the effect of obesity-related traits on COVID-19 severity and hospitalization.

Association of Kidney Comorbidities and Acute Kidney Failure With Unfavorable Outcomes After COVID-19 in Individuals With the Sickle Cell Trait.

Importance: Sickle cell trait (SCT), defined as the presence of 1 hemoglobin beta sickle allele (rs334-T) and 1 normal beta allele, is prevalent in millions of people in the US, particularly in individuals of African and Hispanic ancestry. However, the association of SCT with COVID-19 is unclear. Objective: To assess the association of SCT with the prepandemic health conditions in participants of the Million Veteran Program (MVP) and to assess the severity and sequelae of COVID-19. Design, Setting, and Participants: COVID-19 clinical data include 2729 persons with SCT, of whom 353 had COVID-19, and 129 848 SCT-negative individuals, of whom 13 488 had COVID-19. Associations between SCT and COVID-19 outcomes were examined using firth regression. Analyses were performed by ancestry and adjusted for sex, age, age squared, and ancestral principal components to account for population stratification. Data for the study were collected between March 2020 and February 2021. Exposures: The hemoglobin beta S (HbS) allele (rs334-T). Main Outcomes and Measures: This study evaluated 4 COVID-19 outcomes derived from the World Health Organization severity scale and phenotypes derived from International Classification of Diseases codes in the electronic health records. Results: Of the 132 577 MVP participants with COVID-19 data, mean (SD) age at the index date was 64.8 (13.1) years. Sickle cell trait was present in 7.8% of individuals of African ancestry and associated with a history of chronic kidney disease, diabetic kidney disease, hypertensive kidney disease, pulmonary embolism, and cerebrovascular disease. Among the 4 clinical outcomes of COVID-19, SCT was associated with an increased COVID-19 mortality in individuals of African ancestry (n = 3749; odds ratio, 1.77; 95% CI, 1.13 to 2.77; P = .01). In the 60 days following COVID-19, SCT was associated with an increased incidence of acute kidney failure. A counterfactual mediation framework estimated that on average, 20.7% (95% CI, -3.8% to 56.0%) of the total effect of SCT on COVID-19 fatalities was due to acute kidney failure. Conclusions and Relevance: In this genetic association study, SCT was associated with preexisting kidney comorbidities, increased COVID-19 mortality, and kidney morbidity.

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program.

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (nrs495828 = 53 and nrs505922 = 59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p = 1.32 x 10-199), and thrombosis ORrs505922 1.33, p = 2.2 x10-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10-23, lupus OR 0.84, p = 3.97 x 10-06. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10-13 among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

Host Genetic Liability for Severe COVID-19 Associates with Alcohol Drinking Behavior and Diabetic Outcomes in Participants of European Descent.

Risk factors and long-term consequences of COVID-19 infection are unclear but can be investigated with large-scale genomic data. To distinguish correlation from causation, we performed in-silico analyses of three COVID-19 outcomes (N > 1,000,000). We show genetic correlation and putative causality with depressive symptoms, metformin use (genetic causality proportion (gcp) with severe respiratory COVID-19 = 0.576, p = 1.07 × 10-5 and hospitalized COVID-19 = 0.713, p = 0.003), and alcohol drinking status (gcp with severe respiratory COVID-19 = 0.633, p = 7.04 × 10-5 and hospitalized COVID-19 = 0.848, p = 4.13 × 10-13). COVID-19 risk loci associated with several hematologic biomarkers. Comprehensive findings inform genetic contributions to COVID-19 epidemiology, molecular mechanisms, and risk factors and potential long-term health effects of severe response to infection.

ACE2 Netlas: In silico Functional Characterization and Drug-Gene Interactions of ACE2 Gene Network to Understand Its Potential Involvement in COVID-19 Susceptibility.

Angiotensin-converting enzyme-2 (ACE2) receptor has been identified as the key adhesion molecule for the transmission of the SARS-CoV-2. However, there is no evidence that human genetic variation in ACE2 is singularly responsible for COVID-19 susceptibility. Therefore, we performed an integrative multi-level characterization of genes that interact with ACE2 (ACE2-gene network) for their statistically enriched biological properties in the context of COVID-19. The phenome-wide association of 51 genes including ACE2 with 4,756 traits categorized into 26 phenotype categories, showed enrichment of immunological, respiratory, environmental, skeletal, dermatological, and metabolic domains (p < 4e-4). Transcriptomic regulation of ACE2-gene network was enriched for tissue-specificity in kidney, small intestine, and colon (p < 4.7e-4). Leveraging the drug-gene interaction database we identified 47 drugs, including dexamethasone and spironolactone, among others. Considering genetic variants within ± 10 kb of ACE2-network genes we identified miRNAs whose binding sites may be altered as a consequence of genetic variation. The identified miRNAs revealed statistical over-representation of inflammation, aging, diabetes, and heart conditions. The genetic variant associations in RORA, SLC12A6, and SLC6A19 genes were observed in genome-wide association study (GWAS) of COVID-19 susceptibility. We also report the GWAS-identified variant in 3p21.31 locus, serves as trans-QTL for RORA and RORC genes. Overall, functional characterization of ACE2-gene network highlights several potential mechanisms in COVID-19 susceptibility. The data can also be accessed at https://gpwhiz.github.io/ACE2Netlas/.

Integrative genomic analyses identify susceptibility genes underlying COVID-19 hospitalization.

Despite rapid progress in characterizing the role of host genetics in SARS-Cov-2 infection, there is limited understanding of genes and pathways that contribute to COVID-19. Here, we integrate a genome-wide association study of COVID-19 hospitalization (7,885 cases and 961,804 controls from COVID-19 Host Genetics Initiative) with mRNA expression, splicing, and protein levels (n = 18,502). We identify 27 genes related to inflammation and coagulation pathways whose genetically predicted expression was associated with COVID-19 hospitalization. We functionally characterize the 27 genes using phenome- and laboratory-wide association scans in Vanderbilt Biobank (n = 85,460) and identified coagulation-related clinical symptoms, immunologic, and blood-cell-related biomarkers. We replicate these findings across trans-ethnic studies and observed consistent effects in individuals of diverse ancestral backgrounds in Vanderbilt Biobank, pan-UK Biobank, and Biobank Japan. Our study highlights and reconfirms putative causal genes impacting COVID-19 severity and symptomology through the host inflammatory response.

ACE2 Netlas: In-silico functional characterization and drug-gene interactions of ACE2 gene network to understand its potential involvement in COVID-19 susceptibility.

Angiotensin-converting enzyme-2 ( ACE2 ) receptor has been identified as the key adhesion molecule for the transmission of the SARS-CoV-2. However, there is no evidence that human genetic variation in ACE2 is singularly responsible for COVID-19 susceptibility. Therefore, we performed a multi-level characterization of genes that interact with ACE2 (ACE2-gene network) for their over-represented biological properties in the context of COVID-19. The phenome-wide association of 51 genes including ACE2 with 4,756 traits categorized into 26 phenotype categories, showed enrichment of immunological, respiratory, environmental, skeletal, dermatological, and metabolic domains (p<4e-4). Transcriptomic regulation of ACE2-gene network was enriched for tissue-specificity in kidney, small intestine, and colon (p<4.7e-4). Leveraging the drug-gene interaction database we identified 47 drugs, including dexamethasone and spironolactone, among others. Considering genetic variants within ± 10 kb of ACE2-network genes we characterized functional consequences (among others) using miRNA binding-site targets. MiRNAs affected by ACE2-network variants revealed statistical over-representation of inflammation, aging, diabetes, and heart conditions. With respect to variants mapped to the ACE2-network, we observed COVID-19 related associations in RORA, SLC12A6 and SLC6A19 genes. Overall, functional characterization of ACE2-gene network highlights several potential mechanisms in COVID-19 susceptibility. The data can also be accessed at https://gpwhiz.github.io/ACE2Netlas/.